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Abstract: We investigate cosmological evolution in models where the effective potential 
V ((f)) may become negative for some values of the field <j). Phase portraits of such theories in 
space of variables (4>, (j), H) have several qualitatively new features as compared with phase 
portraits in the theories with V((f)) > 0. Cosmological evolution in models with potentials 
with a "stable" minimum at V(4>) < is similar in some respects to the evolution in models 
with potentials unbounded from below. Instead of reaching an AdS regime dominated 
by the negative vacuum energy, the universe reaches a turning point where its energy 
density vanishes, and then it contracts to a singularity with properties that are practically 
independent of V(<j))- We apply our methods to investigation of the recently proposed 
cyclic universe scenario. We show that in addition to the singularity problem there are 
other problems that need to be resolved in order to realize a cyclic regime in this scenario. 
We propose several modifications of this scenario and conclude that the best way to improve 
it is to add a usual stage of inflation after the singularity and use that inflationary stage 
to generate perturbations in the standard way. 
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1. Introduction 



Since the invention of inflationary cosmology |l]]-[||], the theory of the evolution of scalar 
fields in an expanding universe has been investigated quite extensively, both at the classical 
and the quantum level. While many features of scalar field cosmology are well understood, 
the overall picture remains somewhat incomplete. In this paper we will extend the inves- 
tigation of scalar field cosmology to models with negative effective potentials. We are also 
going to bring together several other issues, such as the impact of radiation and particle 
production on the onset of inflation. This will allow us to get a better understanding of 
various possibilities that may appear in scalar field cosmology. 

We are going to use a general approach based on the investigation of 3d phase portraits 
that show the behavior of the scalar field 0, its velocity and the Hubble constant H = 
We will see that the phase portraits of models with V((/)) > and with V(cj)) < are 
qualitatively different and that additional changes appear when one adds matter and/or 
radiation. 

There are several reasons to study cosmology with negative potentials. The first one is 
related to the cosmological constant problem. The simplest potential used in inflationary 
cosmology is V(4>) = ^m 2 cf) 2 ||]. One can add to this potential a small cosmological 
constant Vq without changing any features of inflation. A small positive Vo ~ 10 -120 (in 
Planck units) would be sufficient to describe the present acceleration of the universe in a 
de Sitter-like state. But why should Vo be so small and positive? What would happen for 
Vo < 0? Does the post-inflationary universe with Vq < behave like anti-de Sitter space, 
which is so popular in M-theory? 

Rather unexpectedly, the answer to this question appears to be negative: After a long 
stage of inflation the universe with Vo < cannot approach an AdS regime; instead of that 
it collapses || [7|, ||. In this paper we will study cosmological behavior in a large class 
of theories with negative potentials and explain why the universe in these theories stops 
expanding and eventually collapses. 

Another reason to study theories with negative potentials is provided by the investiga- 
tion of cosmology in N=2,4,8 gauged supergravity. Recently it was found that in all known 
versions of these theories potentials with extrema at V(<j>) > are unbounded from below. 
Despite this fact, such models can, under certain conditions, describe the present stage of 
acceleration of the universe |?|, || . 

One more reason is related to a formal connection with warp factor /bulk scalar dynam- 
ics in brane cosmology. It has recently been shown that the equations for the warp factor 
and scalar field in brane cosmology with a scalar field potential V(4>) are similar to the 
equations for the scale factor and scalar field in 4D cosmology with the opposite potential 
~V(4>) @- This reveals an interesting relation of cosmology with negative potentials and 
warp geometry with positive potentials. 
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Finally, cosmology with a negative potential V{<j)) is the basis of the recently proposed 



"cyclic universe" model [fR]] based in part on the ekpyrotic scenario [11]. However, unlike 
in the ekpyrotic scenario |Tl]| , the authors of JIO] assume, in accordance with [12], that 
the scalar field potential V(4>) at large (ft is positive and nearly constant. As a result, the 
universe experiences "superluminal expansion" (inflation) that helps to solve some of the 
cosmological problems. In this sense cyclic scenario, unlike the ekpyrotic scenario of Ref. 
[11], is a specific version of inflationary theory rather than an alternative to inflation. 1 Then 



the scalar field rolls to a minimum of its effective potential with V{(p) < 0, the universe 
contracts to a singularity, re-emerges and again enters a stage of inflation. This scenario 



inherits many unsolved problems of the ekpyrotic model [12], including the singularity 



problem [13]. The authors assume that the universe can pass through the singularity and 
that one can use perturbation theory and specific matching conditions at the singularity to 
calculate density perturbations in the post-big bang universe generated by processes prior 



to the singularity [14]. This issue is rather controversial [15|. The possibility of achieving 
a cyclic regime depends on various assumptions concerning the creation of matter and the 
acceleration of the scalar field during the big bang. 

The idea that the big bang is not the beginning of the universe but a point of a 
phase transition is quite interesting, see e.g. ||16||-||22|j. However, the more assumptions 
about the singularity one needs to make, the less trustworthy are the conclusions. In this 
respect, inflationary theory provides us with a unique possibility to construct a theory 
largely independent of any assumptions about the initial singularity. According to this 
theory, the structure of the observable part of the universe is determined by processes at 
the last stages of inflation, at densities much smaller than the Planck density. As a result, 
observational data practically do not depend on the unknown initial conditions in the early 
universe. 

Since the cyclic scenario does require repeated periods of inflation anyway, it would 
be nice to avoid the vulnerability of this scenario with respect to the unknown physics at 
the singularity by placing the stage of inflation before the stage of large scale structure 
formation rather than after it. 

In order to achieve this goal we will examine the conditions that are necessary for the 
existence of the cyclic regime in the model of Ref. j[(J and then check whether the model 
can be modified in a way that would not require various speculations about the behavior 
of matter, the scalar field, and density perturbations near the singularity. 

Our paper will thus consist of two parts. The first part will contain a general study 
of scalar field cosmology with positive and negative potentials. The second part will be 
devoted to a more speculative subject, it will include application of our general results to 
the cyclic scenario. 

In Section [2] we will describe several basic regimes that are possible in scalar field 



1 One should note, however, that this is a very specific kind of inflation that is possible only if the universe 
is exponentially large all the time. Thus the large size of the universe is not explained by inflation in this 
model but rather required for it. 
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cosmology: the universe can be dominated by potential energy, by kinetic energy, by the 
energy density of an oscillating scalar field, or by matter or radiation. The discussion of 
these four distinct regimes will help us to understand the phase portraits of the universe 
that we are going to draw in the subsequent sections. 

Section ^ will describe the use of phase portraits for studying cosmological evolution. 
We will write the evolution equations for the field and scale factor in the form of three first 
order equations plus one time dependent constraint. The solutions to these equations can 
then be represented as trajectories in phase space, clearly showing the possible ways the 
universe can evolve in different situations. Finally, by using a Poincare projection we can 
map the entire phase space onto a finite sphere, thus allowing the complete set of possible 
trajectories to be easily seen. 

In Section ||] we will apply these methods to models with positive definite potentials. 
Such potentials have been extensively studied before with the use of phase portraits ]23|, p4|. 
We study them here partly to introduce the methods we are using and to provide a point 
of comparison for the negative potentials of the following section. We also present some 
new results concerning the effects of matter and radiation on the development of inflation. 

In Section || we show the phase portraits for a model where the effective potential can 
become negative. We discuss general properties of such models, and in particular the ways 
in which they differ from the models of the previous section. One of our major conclusions 
is that such models generically enter a stage of contraction. In Section [6| we will examine 
in detail the transition from expansion to contraction in models of this type. 

Many of the features of scalar field cosmology that we are going to discuss are model- 
independent. The phase portraits in Sections ||-|| all use the simplest model V{<j)) = 
m 2 (j) 2 /2 + Vb, but in Section ^ we discuss some other theories with negative potentials. 

In Section |8| we will discuss cosmological evolution near the initial and final singulari- 
ties, and in particular the role of particle production and anisotropy near the singularity. 

In Section ^ we will apply our methods to the investigation of the cyclic scenario. As 
we will see, the cyclic regime in this scenario does not appear automatically. One should 
fine-tune the potential V(4>) and learn how to work with the super-Planckian potentials 
|V(0)| > 1. One should also introduce superheavy particles with specific properties, study 
their production at the singularity, and make sure that they do not affect the present stage 
of the evolution of the universe. This adds new "epicycles" to this scenario, making it even 
more speculative. We discuss several possible modifications of this scenario and conclude 
that the best way to improve it is to add a usual stage of inflation before the stage of large 
scale structure formation. This modification resolves many problems of the original version 
of the cyclic scenario. In this modified form of the cyclic scenario inflation is once again 
the source of density perturbations as well as the resolution of the cosmological problems 
such as homogeneity and flatness. 

Section |l^ summarizes our main conclusions concerning cosmology with negative po- 
tentials and cyclic universe. 
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2. Four Basic Regimes in Scalar Field Cosmology 



2.1 A toy model with V{<j>) = \m 2 4> 2 + V 

We will study the behavior of a homogeneous scalar field in a Priedmann universe with the 
metric 

ds 2 = -dt 2 + a 2 {t)dsl , (2.1) 
where ds^ = r )ijdx l dx 3 is the metric of a 3d space with constant curvature, k = 0, ±1. 

In this paper we will use a system of units in which M p = 1, where M p = (&irG)~ l l 2 ~ 
2 x 10 18 GeV. The Friedmann equation for a scalar field with potential energy density V(<f>) 
is 

*" = (;)' = 5'-? = 5 G* +v( * )+fc H- (2 - 2) 

Here p is the total energy density and p a is the density of matter with equation of state 
p a = ap a . For non-relativistic matter a = 0, while for radiation a = 1/3. 

The evolution of H is given by a combination of the Einstein equations 



H = -\(p + p) + \ = -\{^ + p a (l + a)) + ^ . (2.3) 



Alternatively, one can use the equation 

- = ~\(P + 3P) = \(V{4>) - 4?) - \pa{l + 3a) . (2.4) 
a b 6 b 

The evolution of the scalar field <j) follows from the Einstein equations, 

4> + 3H<j) + V <t> = . (2.5) 

We shall study the basic properties of scalar field cosmology using as an example the 
simplest harmonic oscillator potential 

V{ij>) = \m 2 ^ 2 + V , (2.6) 



Surprisingly, we will find that cosmology with the potential ( |2.q ) with Vo < shares some 
common features with the cosmology of the "inverse" harmonic oscillator potential 

V{4>) = ~\m 2 4> 2 ~ Vo . (2.7) 

In particular, the expansion of the universe in theories with Vo < always turns into 
cosmological contraction. 

Constructing phase portraits is a powerful method for investigating the dynamics of 
the scale factor/scalar field system (|2 . 3| ) - (|2 . 5|) . Before we look at the phase portraits for 
various values of Vo in this model, it will be useful to discuss some of their features. For 
the remainder of this section we will consider k = 0, i.e. flat universes. While this case will 
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be the main focus of our discussion throughout the paper, we will in several cases refer to 
the extension of our results to open or closed universes as well. 

There are four basic regimes that we may encounter: the universe can be dominated 
by the potential energy density V(<f>), by the kinetic energy density (f) 2 /2, by the energy 
density of an oscillating scalar field, in which case V(cfr) ~ <ft 2 /2, or by matter/radiation 

Pa- 



2.2 The inflationary regime: Energy density dominated by V{4>) 

Inflation occurs when the energy density is dominated by V{<f>). In this case <p 2 /2,p a -C 
V{4>) an d |^| <C |3.ff</>|. This corresponds to the vacuum-like equation of state 

p = -p. (2.8) 

The equations for a and <p in this regime have the following form: 

2 /d\ 2 mV V . . 

H =UJ = -f + Y- (2 - 9) 

3-<j) + m 2 (j) = . (2.10) 
a 

The solutions of the equations for <j)(t) and a(t) for the most interesting case — ' £ S> \Vo\ 
are given by 0. [2^] 

<Kt) =<t>o~ \[\rnt , (2.11) 

a(t) = a exp . (2.12) 

These solutions, which describe inflationary expansion, are valid only for <fi 2 /2 <C V((j>), 
which implies that inflation ends at 

\(f>e\ ~ 1 • (2.13) 

2 i 2 

In this paper we will assume that m 2 3> \Vo\, in which case m 2 3> \Vq\ is always satisfied 
during inflation. 

Note that the same solution is valid if one reverses the time arrow, t — > —t, in which 
case it describes a quasi-exponential contraction of the universe (deflation). 



2.3 The kinetic regime: Energy density dominated by (f> 2 /2 

Another important regime occurs when the energy density is dominated by <j) 2 /2. In this 
case V(4>),p a -C 4> 2 /2 and |<ji>|, |3if(/>| 3> m 2 4>. This corresponds to the "stiff" equation of 
state 

p = p. (2.14) 
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The equations for a and 4> are: 

(ha , 

X = -3- . (2.16) 

cf> a 

The solutions can be written as follows: 

a{t) = t 1/3 , (2.17) 

*-*>=■=)/! >«t ! f-^- (2 ' 18) 

These solutions can describe an expanding universe or a universe collapsing towards a 
singularity. 

During the expansion of the universe, the inflationary regime V((p) 3> </> 2 /2 represents 
a stable intermediate asymptotic attractor. Even if a flat universe begins in a state with 
(j) 2 /2 3> V(4>), it typically rapidly switches to an inflationary regime with V((p) 3> <fi 2 /2 
1 23 , |24], 53], This occurs because during the expansion of the universe with </> 2 /2 3> ^(</>), 
the value of the kinetic energy drops down like t~ 2 , whereas the field changes only loga- 
rithmically. Therefore for all power-law potentials, the value of V{4>) decreases much more 
slowly than 4> 2 /2. When it becomes greater than c/> 2 /2, inflation begins. 

During the collapse of the universe, the opposite occurs. V{4>) grows only logarithmi- 
cally, whereas <^> 2 /2 diverges as t~ 2 , where t is the time remaining before the big crunch 
singularity. This means that the regime <ft 2 /2 3> V((f>) generically occurs at the stage of 
collapse. In this regime one can neglect V((f>) in the investigation of the singularity at 
t -» 0. 



2.4 The oscillatory regime: Evolution determined by the energy density of an 
oscillating scalar field 

Now let us assume that the field <p oscillates near eft = with frequency much greater 
than H, and that the average value of V(<j)) during these oscillations is much greater than 
Vo = V(0). In this case one can neglect the term 3Hcj) in Eq. ( |2.5| ), so that in the first 
approximation one simply has 

4> + m 2 <p = (2.19) 

and 

4> = $ s'mmt . (2.20) 

Here $ is the amplitude of the oscillation. The pressure p = (p 2 /2 — V(4>) produced by these 
oscillations is given by ^$ 2 cos 2mt, so if one takes an average over many oscillations, 
the pressure vanishes, p ~ 0. The universe in this regime expands as a ~ Since the 
total energy of pressureless matter is conserved, the amplitude of the oscillations decreases, 
~ a" 3 / 2 ~ r 1 . 
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The regime of oscillations usually begins after the end of inflation, at <j> ^ 1. As long as 
one can neglect Vq, the field oscillations after inflation approach the following asymptotic 
regime p7||: 

2\/2 . V2 



sin mt 



sin mt . 



y/Smt N 
Here t is the time after the end of inflation and N is the number of oscillations. 



(2.21) 







a 4i 



lOMf JDOO 3000 4000 5000 



1000 2000 3000 4000 5000 



Figure 1: Evolution of the scalar field and the scale factor in the model V(4>) = ^-0 2 + Vb with 
Vo > 0. In the beginning we have a stage of inflation with the field <fi linearly decreasing at <f> > 1- 
At this stage the equation of state is p w —p. Then the field enters a stage of oscillations with a 
gradually decreasing amplitude of the field; p -C p. When the energy of the oscillations becomes 
smaller than Vq, the universe enters a second stage of inflation. 

It is amazing that this simple model with Vo > can describe not only inflation in 
the early universe, but also the present stage of inflation/acceleration. Indeed, when the 
amplitude becomes very small the term Vo will become important, and the universe enters 
a second stage of inflation with H 2 = Vq/3. The amplitude of oscillations of the field <f) in 
this regime falls down exponentially. In particular, for m 2 ^> H 2 the amplitude decreases 
as e~ 3Ht / 2 . The evolution of the scalar field and the scale factor in the theory with Vo > 
is shown in Fig. |l|. 

Meanwhile, if one considers the model with Vo < 0, a dramatic change occurs when the 
energy density of oscillations (and matter) gradually decreases and becomes comparable to 
— Vq. According to Eqs. (|2.2[) and (|2.3|) , the expansion of the universe slows down at that 
time, and eventually the universe begins collapsing, see Fig. ||. 

When the universe contracts, the amplitude of oscillations grows as a~ 3 / 2 . However, 
this process does not continue too long. Indeed, let us compare 3H(p and m 2 (f) in this 
regime. If one can neglect Vq (and this is always the case for a sufficiently large one 
has H ~ m$/\/6 and <fi ~ m<j). Therefore one has |3iJ(^| 3> |m 2 0| for <p S> 1, so instead of 
Eq. ( 2.2C| ) one should use Eq. ( p. 16 ). Thus, during the collapse of the universe the stage 



of oscillations ends and the regime dominated by kinetic energy begins at 

\<k\ ~ 1 • (2-22) 

Note that \(j) b \ ~ \(j) e \, see Eq. (pl^ ). 
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Figure 2: Evolution of the scalar field and the scale factor in the model V{4>) = ^-0 2 + Vb with 
Vo < 0. In the beginning we have a stage of inflation with the field <fi linearly decreasing at <f> > 1- 
At this stage the equation of state is p w —p. Then the field enters a stage of oscillations with a 
gradually decreasing amplitude of the field; p -C p- When the energy of the oscillations becomes 
equal to | Vb|, the universe stops expanding and begins to contract. At this stage the amplitude 
of oscillations grows. When it becomes greater than 0(1), the field stops oscillating, the energy 
density is dominated by the kinetic energy of the scalar field, p ~ p, and the universe collapses. 



We will study the switch from expansion to contraction in a flat universe in a much 
more detailed way in Section ||. However, we would like to make here some comments 
concerning this process. 

The general textbook wisdom is that open and fiat universes expand forever, whereas 
closed universes eventually collapse. This lore was based on investigation of universes 
with vanishing cosmological constants. A closed universe with a sufficiently large positive 
cosmological constant may expand forever, whereas open and flat universes with a negative 
cosmological constant eventually collapse. 

One of the well-known solutions of this type is an open universe with a negative 
vacuum energy Vq. There is a solution to the Friedmann equation H 2 — a~ 2 = Vn/3 for 

Vq < 0: a(t) = Y^j^y sin \J This is a specific section of anti de Sitter space, which 
is popular in M-theory and brane cosmology. This universe has a coordinate singularity 
at t = 7r^/|^|-. Naively, one might think that this is exactly what we have found in our 
investigation of universes with Vo < 0, namely that when the energy density of matter 
in an expanding universe decreases and the total energy density becomes dominated by 
a negative cosmological constant, our universe reaches an AdS regime dominated by a 
negative cosmological constant. 

However, this is not the case. We discuss here a flat universe regime, which appears 
after a long stage of inflation. In this case (unless one considers open inflation models with 
< 1) the term \ with k = ±1,0 can be omitted in the general Friedmann equation. The 
Friedmann equation H 2 = p/3 describing a flat universe does not have any solutions with 
p{4>) < 0. Once the universe approaches the turning point where the total energy density 
vanishes it begins collapsing, and the total energy density becomes positive again [||, [j], 
Thus the standard inflationary prediction f2 = 1 implies that we cannot live in AdS space 
dominated by a negative cosmological constant [|], [|. 
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2.5 Evolution determined by the energy density of matter or radiation 



The first models of inflation were based on the assumption that the universe from the very 
beginning was in a state of thermal equilibrium; inflation began when the temperature of 
the universe became much smaller than the Planck temperature T ~ M p |2], Later it 
was found that this assumption is not necessary, and in many models inflation may start 
immediately after the big bang [Q. In this case the existence of matter prior to inflation 
becomes less important, and sometimes it even hampers the development of inflation p5| . 
Therefore many works on initial conditions for inflation neglect the possible impact of mat- 
ter on the motion of the scalar field and concentrate on finding self-consistent cosmological 
solutions describing scalar fields in otherwise empty universes. This is the simplest ap- 
proach, especially in cases where o 2 /2 <C V{cft) and inflation begins immediately after the 
big bang. 

However, in some cases the scalar field initially may have large kinetic energy, ^ 2 /2> 
V{(ft). Moreover, one may expect creation of relativistic or non-relativistic particle near 
the singularity. Note that the existence of even a small amount of matter may have an 
important effect on the motion of the field. Indeed, the kinetic energy of the scalar field 
<ft 2 /2 in the regime (ft 2 /2 3> V((/>) decreases as a -6 . Meanwhile, the density of radiation 
decreases as a~ 4 and the density of non-relativistic matter decreases as a -3 . Therefore the 
energy density of matter eventually becomes greater than <ft 2 /2. As we will see, once it 
occurs, the field rapidly slows down or even completely freezes. This effect may provide 



good initial conditions for a subsequent stage of inflation [28|. 

Indeed, let us assume that in the beginning the field (ft moves very fast, so that |3J3"0| S> 
|V,0 | = |m 2 (/>|. Suppose, however, that at some moment the energy density of the universe 
becomes dominated by matter with the equation of state p a = ap a . In this regime one can 



represent the cosmological evolution in the following form [25] 

-3(l+ex) 

Pa = Pa{t\ 



a(t) = a 
H 





'a(t)\ 


o)| 






v a / 




2 




, 3(l+a) 


toy 




2 





2 

1 + a 



(2.23) 



3(l + a)t ' 

l X a o X f to 
9 = = <Po — 

a 6 \t 

This regime has a very interesting feature: Even if it continues for an indefinitely long 
time, the change of the field (ft during this time remains quite limited. Indeed, 

OO OO 2 

A<ft < J fat = <ft J (j) ^ dt = <ft t . (2.24) 



to to 

iof inn ( /O , n \ + In on A\^-f-„ ^ , „ 



If to is the very beginning of matter domination (4^/2 ~ p a ), then o io ~ -7^77 — r = 0(1). 
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Therefore 

A0 £ 1 (2.25) 

in Planck units (i.e. A(j) <^ M p ). This means, in particular, that a free field <f> in a matter 
dominated universe cannot move by more than 0(M p ). 

This simple result has important implications. In particular, if the motion of the field 
in a matter-dominated universe begins at \cj>\ 3> 1, then it can move only by A<j) ^ 1. 
Therefore in theories with flat potentials the field always remains frozen at \<p\ 3> 1. 

The field begins moving again only when the Hubble constant decreases and |3.ff0| 
becomes comparable to \V,<p |. But in this case the condition 3H(j) « \V„p | automatically 
leads to inflation in the theory m 2 4> 2 /2 + Vo for \Vq\ < m 2 and ^> 1. 

This means that even a small amount of matter or radiation may increase the chances 
of reaching a stage of inflation, see [2S] and Fig. ||| in Section [|. Indeed, consider any 



theory with V((p) ~ 4> n . Suppose in the beginning we had a kinetic energy dominated 
regime <fi 2 /2 p a ,V(<f>) starting at <fi S> 1. Then the field 4> would change very slowly, 
whereas (j) 2 /2 would rapidly drop down until it became comparable either to V(<p) or to 
p a . If at that moment V(4>) > p a , inflation would begin immediately. But even in the 
most unfavorable case V(cj)) <C p a inflation would begin eventually. Indeed, at 3> 1 one 
has the double inequality m 2 = V" <^ V(4>) <C p a ~ H 2 . Therefore the Hubble constant 
is much greater than the effective scalar field mass. In this case the field practically does 
not move until the desirable regime V(4>) > p a is reached and inflation begins. 



3. Phase Portraits and Cosmological Evolution 

Having discussed some important limiting regimes in scalar field cosmology, we are now 
ready to investigate the complete evolution of a Friedmann universe with a scalar field. 
Later we will discuss the effects of adding matter to this system, but for now we restrict 
ourselves to a system with three independent variables, <p, <p, and H. To study this system 
we find it most convenient to rewrite the evolution equations for a and <j) as a set of three 
coupled, first-order, differential equations: 



dt 



(3.1) 



f = -3^ - V >(t> (3.2) 

dH 1 t 12 tA u2 



plus the constraint equation 



dt ,> V)-H' (3.3) 



H 2 - If -\v= -4r. (3.4) 
6 3 cr 



All solutions to these three equations can be represented as trajectories in the 3d 
phase space of (f>, 4>, and H. Simply looking at plots showing a number of these trajectories 
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can help give some intuition for the cosmology of a particular model (as defined by the 
potential V). There are a number of ways to get more information out of the phase 
portraits, however. 

One important step is to determine all of the critical points, i.e. the points for which 
the derivatives of all three phase variables vanish. There are finite and infinite critical 
points. Every trajectory must begin and end at these critical points. 

To find infinite critical points and visualize the flow of trajectories at infinity, a useful 
trick is to do a Poincare mapping 

where x is any of ((j), <fi, H) and r 2 = (ft 2 + (ft 2 + H 2 . The interior of the unit sphere 
4>p + 9p + Hp = 1 maps to the infinite phase space of (f>, (j), and H, so by plotting 
trajectories in these new coordinates the entire phase space can be easily visualized. At 
times in this paper we will plot a 2d phase portrait, e.g. in the variables 4> and <p only. In 
these cases we use a 2d Poincare mapping where r 2 = (j) 2 + 4> 2 . 

With the Poincare mapping it is possible to identify a set of infinite critical points, 
namely those that occur on the bounding sphere (j)p + (pp + Hp = 1. These points represent 
the possible starting and ending points for all trajectories that go off to infinity in the usual 
coordinates. 

Because no two trajectories can ever cross in phase space, it is easy to define the behav- 
ior of a system whose phase portrait is two dimensional. Fortunately, for the cosmological 
systems we are considering we can identify a 2d surface that separates different regions of 



the 3d phase space. For the flat universe k = the constraint equation (3.4) defines a 2d 
surface. All trajectories in this case are located at this surface, i.e. the phase portrait for 
the flat universe is two dimensional. This surface in turn divides the phase space into three 
separate regions (including the surface itself) representing the possible types of curvature. 
No trajectory can pass from one of these regions to another. Although the location of the 
finite critical points for a given model depends strongly on V, the structure of the infinite 
critical points is very similar across a wide range of potentials. See || for recent discussion. 

4. Cosmology with a Non-Negative Potential 



As a simple example we consider the model V{<j)) = V$ + \m 2 4> 2 discussed in Section 2.1 



By rescaling the field and time variables the mass m can be eliminated from the equations, 
so for simplicity we simply set m = 1 in what follows. Thus the evolution and constraint 
equations become 

^ = -34>H - 4, (4.2) 
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6H 2 -ft -<j) 2 - 2V = -6-s . (4.4) 

or 

The hypersurface representing a flat universe is given by setting k = in the constraint 
equation, which gives 

m 2 -ft -ft = 2V . (4.5) 

The surface defined by this equation is a hyperboloid. For positive definite potentials 
Vo > it is a hyperboloid of two sheets, meaning the two branches at H > and H < 
are disconnected. For Vo = this hyperboloid reduces to a double cone. 

There are two finite critical points for this system at <j> = § = 0, H = ± For 
Vo = these two points reduce to a single finite critical point at the origin. To find the 
infinite critical points we first rewrite the evolution equations in terms of the Poincare 
variables and then set their derivatives equal to zero. This yields eight points. 

Figure ^ shows the phase space for this model with Vo > along with a sample of 
trajectories for k = 0. The hyperboloid along which all of these trajectories lie represents a 
flat universe. The upper branch corresponds to expansion and the lower one to contraction. 
The fact that the two branches are disconnected means that in a flat universe in this 
model expansion can never reverse and become contraction. Note that this conclusion is 
unchanged for the case Vo = 0. In that case the hyperboloid becomes a double cone and 
the two branches touch at a single point. Since that point is a critical point, however, 
no trajectories can pass from one branch of the cone to the other. The lower branch 
corresponds to the upper branch with time reversal t — * —t. The upper branch of the 
flat universe hyperboloid is shown projected into a 2d plot in Figure [|. 2 This plot is very 
similar to the one shown in |23[] for this model with Vb = 0. 

For an expanding universe there are four infinite critical points, two repulsors labeled 
Ri and R2 and two saddle points labeled Si and S2. All trajectories begin at R\, R2 and 
wind towards the focus at the center. The separatrices emanating from S% and S2 represent 
attractor trajectories (not to be confused with attractor critical points). Along these tra- 
jectories the universe experiences inflation (ft 3> ft) until it nears the center and begins 
winding around it, corresponding to field oscillations near the potential minimum. These 
separatrices represent a set of measure zero in the space of trajectories; the two shown are 
the only trajectories that begin at the saddle points. Nonetheless they are important be- 
cause most of the trajectories emanating from the repulsor points asymptotically approach 
the separatrices. This is why inflation is a generic feature of models such as this one, and 
also why inflation erases all information about the initial conditions that preceded it. 

Thus a typical trajectory passes through three of the four regimes described in Section 
|||. Near the repulsors the kinetic energy dominates and the equation of state is stiff, p p. 
Near the main part of the separatrices the equation of state is inflationary, p ~ —p. Finally 
near the center the scalar field oscillates and the equation of state is that of non-relativistic 

2 This is not a direct "shadow" since it uses the 2d rather than the 3d Poincare mapping, see Section 
Effectively the upper branch of the hyperboloid is stretched out onto the circle rather than vertically 
projected down to it. From here on we will refer to such 2d portraits as projections of the 3d ones. 
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Figure 3: Phase portrait for the theory V((f>) — ^m 2 <j) 2 + Vq with \q > in rescaled coordinates 
(0, </>,-ff). The branches describing stages of expansion and contraction (upper and lower parts of 
the hyperboloid) are disconnected. 



matter, p <C p- During the oscillations the scalar field decreases as 

4>(t) ~ j^: sin mt , (4.6) 

where N is the number of oscillations, see Eq. ( p. 21 ). Although particle production is 
not included in these phase portraits, this evolution will typically end with the scalar field 
decaying into other forms of matter, thus finishing the evolution in the fourth regime, 
matter and/or radiation domination. The contracting branch is a mirror image of the 
expanding one, with the same three regimes occurring in the opposite order, finally ending 
with a big crunch singularity at the attractor points Ai and Ai- 

For an open or closed universe the trajectories would lie in the interior or exterior 
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Expanding Universe 
Rl 




Figure 4: Projection of the upper branch of the full phase portrait for the theory V(4>) = im 2 </> 2 
V with V > in rescaled coordinates {(f), <p). 



of the hyperboloid, respectively [23|. For an open universe nearly all trajectories would 

asymptotically approach the separatrices on the flat universe hypersurface. This tendency 

reflects the fact that for most initial conditions inflation will occur and drive the universe 

towards flatness. Once this has occurred the trajectories spiral in towards the focus at 

the bottom of the hyperboloid. For a closed universe there are also many trajectories that 

rapidly approach these separatrices, but there is also a class of trajectories that moves 

from the repulsive critical points to the attractive ones without ever passing near the 

flat universe hypersurface. These trajectories reflect closed universes that collapse rapidly 

before inflation has a chance to occur. 

This conclusion becomes even more apparent if one takes into account matter /radiation 

|28[| . As we have argued in Section 2Jj , the existence of matter rapidly freezes the motion 

of the scalar field. Therefore if the field <f> was initially large and had a large velocity such 

that (f) > 1, 4> 2 /2 > V(<j>), then the presence of matter would increase the probability of 

inflation. This can be confirmed by comparing the phase portraits of the universe with and 

without radiation. Although the phase portrait with radiation is three dimensional, it is 

convenient to make its projection to the (<p, cj>) plane, see Fig. [5|. 

In the second and fourth quadrants of this figure the field starts out moving towards 
the minimum. The presence of radiation slows the field down, causing it to move more 
quickly towards the inflationary separatrix trajectory. In the first and third quadrants 
where the field starts out moving away from the minimum the duration of inflation is 
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4> 




Figure 5: Phase portrait for the theory V(4>) = ^m 2 4> 2 without Poincare mapping. The blue 
(thick) lines show trajectories describing the universe without radiation. The scalar field has half 
Planck density at the beginning of the simulations. The red (thin) lines show trajectories where an 
equal amount of energy in radiation was added to the system. As we see, in the presence of radiation 
the velocity of the scalar field rapidly decreases, which usually leads to the onset of inflation. 



slightly diminished by the presence of radiation, but the probability of inflation is nearly 
unity. 



5. Cosmology with a Negative Potential 

Now we turn to the main subject of our investigation, cosmological models with scalar 
field potentials that may become negative. We will continue using the simple example 
V(<f>) = V + \m 2 (t) 2 , but now we will consider Vq < 0. The hyper surf ace representing a 
flat universe is still defined by 

6H 2 -<]) 2 -0 2 = 2V , (5.1) 

but with Vq negative, this surface is a hyperboloid of one sheet. 

Figure ^ shows the phase space for this model and sample trajectories for a flat universe. 
The phase space is two dimensional, but its topology is very different from that for non- 
negative potentials. The infinite critical points are unchanged because the finite term Vq 
has no effect at infinity, but there are no finite critical points. Thus all trajectories begin 
at infinity with H > and end at infinity with H < 0. This is possible because the regions 



corresponding to expansion and contraction are now connected. This property is valid for 
all types of curvature k, i.e. for open, flat or closed universes. 




Figure 6: Phase portrait for the theory V(</>) = \m?(f) 2 + Vo for Vq < 0. The branches describing 
stages of expansion and contraction (upper and lower parts of the hyperboloid) are connected by a 
throat. 

To show a 2d projection of the flat universe hypersurface for this model, we have to 
plot both the expanding and contracting branches, as depicted on Figure ^. Trajectories 
in the expanding universe region spiral in towards the center. When they touch the inner 
circle, the "throat" of the hyperboloid, they pass into the contracting universe region. 
There they spiral back out to infinity, i.e. the big crunch. Thus typical trajectories in 
this scenario pass through the three regimes described above, kinetic energy domination, 
potential energy domination, and oscillations, and then pass back through them in reverse 
order. As before, including particle production will typically introduce a matter /radiation 
dominated regime after the first stage of oscillations. Eventually, however, the matter and 
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radiation will redshift away and the universe will begin contracting. We will examine this 
process in more detail in the next section. 




Aside from this "wormhole" connecting the expanding and contracting branches this 
phase portrait looks a lot like the one for Vo > shown in Figure ||. Note, however, that 
in this case the separatrices emanating from the saddle points S\ and S2 no longer spiral 
in to the center, but rather end up reaching the points A\ and A2. Likewise there are 
separatrices that begin at R\ and R2 and end on S3 and S4. In the expanding phase their 
segments and segments of nearby trajectories represent the rare cases that manage to avoid 
inflation. In the contracting phase they become the marginal trajectories separating those 
that end at positive and negative <j>. The number of windings (i.e. field oscillations) can 
be estimated by setting m 2 (j) 2 /2 = \Vq\ and using ( |4.6[ ) to give 

(5.2) 




(This number of windings can be used to determine which repulsors and attractors are 
connected to which saddle points, e.g. whether the separatrix that begins at Ri ends at 
S3 or S4.) 

The phase portraits shown above were constructed in a way symmetric with respect 
to time reversal, t — > —t. This is a legitimate approach, since our equations allow all 
of the solutions shown in the previous figures. However, one can obtain some additional 
information if, for example, one considers trajectories equally distributed with respect to 
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Figure 8: A different version of the phase portrait for the theory V((f>) — ^m 2 ^ 2 + Vo for Vq < 0. 
We begin with the trajectories evenly distributed with respect to the initial values of 4> in the early 
universe (upper part of the hyperboloid) and see what happens to them in the lower part. These 
trajectories are concentrated near the red separatrices and repulsed from the green ones. 

the initial value of the field (ft at the Planck time and follows their evolution from the region 
with H > to the region with H < 0. 

If we do so, the phase portrait shown in Fig. |6] starts looking somewhat different. 
Almost no trajectories beginning in the upper part of the hyperboloid are seen in its lower 
part, and those few that can be seen there are positioned very close to the (red) separatrices 
going from 5*1 to A2, and from 52 to Al see Fig. ||. No trajectories are seen near the 
(green) lines going from Rl to 53 and from R2 to 54. This might seem surprising because 
these lines are solutions of the equations of motion, so there must be other solutions nearby. 
Indeed we have seen them in Fig. [?]. However, the (red) lines going from 51 to A2 and from 
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52 to Al are strong attractors in the regime H < 0, whereas the lines going from Rl to 

53 and from R2 to 54 are strong repulsors. Therefore most of the trajectories originating 
at H > and homogeneously distributed with respect to the field (ft a t the Planck density 
are repelled from the lines going from Rl to 53 and from R2 to 54, and tend to merge 
with the lines going from 51 to A2 and from 52 to Al. 

This effect is especially apparent in the 2d phase portrait, where we do not make the 
Poincare mapping, see Fig. ||. Most of the trajectories coming from the panel with H > 
have merged with the red separatrix on the panel corresponding to H < 0. 



Expanding Universe 



Contracting Universe 




4> 




4> 



4> 



Figure 9: As in the previous figure, we begin with the trajectories evenly distributed with respect 
to the initial values of <p in the early universe. However, now we show a 2d projection of these 
trajectories, without Poincare mapping. Left: {<p, (f>) projection of the H > branch. Right: (0, <fi) 
projection of the H < branch. Trajectories from the left panel continue on the right panel. 

An important (and obvious) feature of the 3d phase portraits Fig. ^ and Fig. || 
is that the separatrices, as well as other trajectories, never intersect in 3d. This is a 
trivial consequence of the fact that we are solving a system of 3 first order equations for 
3 variables, 4>, <j> and H. One of implications of this fact is that a bunch of trajectories in 
the immediate vicinity of the (green) lines going from Rl to 53 and from R2 to 54 never 
reach the inflationary regime described by the (red) inflationary separatrices going from 
51 to A2 and from 52 to Al. Only the trajectories that are sufficiently far away from the 
green lines going from Rl to 53 and from R2 to 54 can enter the stage of inflation. 



This observation will be important for us when we describe the cyclic scenario ]10[, 
see Section |9| In this regime the red inflationary separatrices reach the singularity and 
are supposed to bounce back. In the language of the phase portraits this bouncing back 
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implies that the end of the red line going from SI to A2 becomes the beginning of the 
green line going from Rl to 53. But in this case the universe cannot attain the inflationary 
regime, since the trajectories close to the green line never switch to the vicinity of the red 
line. Thus the cyclic regime is possible only if bouncing from the singularity shifts the 
trajectory to the right from the green separatrix. From Fig. |9] it is obvious that this shift 
may happen either due to an increase of <fi or due to an increase of the field <f>. 

The evolution of this system in an open or closed universe is not very different from 
the flat universe evolution, although the phase space is three dimensional. Because of 
the structure of the trajectory flow between their ends at the infinite critical points, all 
trajectories pass from expansion to contraction, even for an open universe. As with Vq > 
the trajectories for the open and closed cases will tend to asymptotically approach the 
flat universe hypersurface, and more specifically will tend to approach the inflationary 
separatrices. As before, however, the closed universe will include some trajectories that 
quickly collapse before experiencing inflation. 

It is instructive to estimate the time that the universe may spend in its post-inflationary 
expanding phase before it begins to contract. The energy density of the oscillations of the 
scalar field, just like the energy density of nonrelativistic matter, decreases as pcdm ~ 
The universe begins to collapse at pcDM+^b = 0. This happens at t ~ , 2 -- . As one could 

expect, this time can be greater than the present age of the universe only if |Vo| <^ 10~ 120 . 

This estimate remains true for a wide variety of potentials and for matter with any 
reasonable equation of state. However, in the theories where V(4>) has a very flat plateau 
or a local minimum, the universe may spend a very long time before the field 4> falls down 
to the minimum with V((f>) < [0, [l(J. Therefore in general the life-time of the universe 
may be very large even in theories with a very deep minimum of V{<j)). 

6. Going from expansion to contraction in the model V(<fi) — ^(p 2 + Vq 

Having analysed general properties of phase portraits in the theory V((j)) = ^-<j) 2 + Vq, let 
us study in a more detailed way the most interesting feature of the models with Vo < 0, 
the switch from expansion to contraction. It is always possible to study this process 
numerically, but sometimes one can do better than that. 

It will be convenient to represent V{<p) = ^%-<p 2 + Vq in the form 



V( ( f>) = 




This potential has a minimum at (j) = 0, where it takes a negative value V{4>) = —^<pQ. 
The potential vanishes (V(4>) = 0) at <p = ±00 • 

Let us assume, in the first approximation, that the scale factor of the universe does 
not change much during each oscillation of the field (p. In such a case the field <ft would 
experience a simple oscillatory motion, 

(f)(t) = $ cosmt , (6.2) 
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where <3? is the amplitude of the oscillations. In this case the total energy density of the 
scalar field would remain constant, p = ^(<£ 2 — </>q). 

This approximation works well for $ ~ <pQ. For <I> > 0o, there are two cosmologi- 
cal solutions, describing either an expanding universe with H = +m v / ($ 2 - $j)/6 or a 
contracting universe with H = -mi/ (<£ 2 — <^>o)/6. 

If the Hubble constant H is positive, the amplitude of the field and its total energy 
density decrease. If the initial amplitude of the oscillations is much greater than <fo, the 
field oscillates with a slowly decreasing amplitude until it approaches 4>o. But the energy 
density cannot decrease too much because at the moment when p = V((p) + (p 2 /2 vanishes, 
the Hubble constant vanishes too, so that d = 0. Then the universe begins to collapse, 
d < 0, and the amplitude of the oscillations begins to grow. Eventually this growth becomes 
so fast that the field stops oscillating and moves towards <p = ±00. 

The best way to understand this effect is to examine what happens during the critical 
oscillation when the sign of d changes. We will study this process analytically, making 
some simplifying approximations. 

First of all, we will assume that the field <f> begins this oscillation at t = moving with 
zero initial velocity from a point <fii « <po such that < A(f> = (pi — (po <C (po- The initial 
energy density of the field is AV = V{(pi) = -£(<£i - <#j) <C |^(0)|. We will try to evaluate 
the turning point moment t c where d = (i.e. H = 0). 

Let us consider the series expansion of the Hubble parameter around the beginning of 
this process 

H(t) *Hi + Hit + ^Hit 2 + ^Hit 3 + ... , (6.3) 

where Hi and its derivatives are taken at t = 0. The reason to include the terms up to 
t 3 in this series is the following. From the relation H = —\(j) 2 we find that for vanishing 
initial velocity tfii = one has Hi = Hi = 0. The first nonvanishing coefficient Hi ~ —(f) 2 ~ 
-{V'{4>i)) 2 = -m A 4>\ is negative. Note that Hi = y/V(4>i)/3 = This means that 

at the moment 



(i"W>i))V V'" 2 *o/ 

the Hubble parameter vanishes. Note that the first part of this equation is pretty general, 
whereas the second one is specific to quadratic potentials. 

At the turning point 

3AV ^ 1/3 



( 2/; >o 



(6.5) 



These results imply that the turn occurs during the first oscillation starting at (pi if AV^ <^ 



rn 



: 1 i.e. (f)\ — (pQ ^ <pQ. In the most interesting case AV^ <C m 2 4>Q the turn occurs in the 



immediate vicinity of the point 4>q where the potential becomes negative. 
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To study the subsequent evolution of fat) and a(t), let us assume that the scale factor 
a during the first oscillation does not change much. This is a reasonable assumption since 
a = at the turning point. We will therefore take a = 1 during this oscillation, and 
fat) = fa cos mt. The potential energy density of the field is 



V(fa = AV- sin 2 mt (6.6) 

and the acceleration of the universe is given by 

.. a V-<j> 2 AV m 2 fa , , w , 

a w - = — = sin 2 mt . 6.7) 

a 3 3 2 v ; 

Taking into account that initially a = a 7AT73 « yAy/3, this yields 

t P-t + sin 2mt + ^AV/3 . (6.8) 



3 4 

By integrating this relation from t = to t = n/m, i.e. during one half of an oscillation, 
one finds that the condition a ~ 1 implies then that </>i ~ fa <C 1, i.e </>o <C M p . 

Now we are going to find how the energy density p of the field (j> changes during the 
time ir/m when the field (j) moves from fa to —fa. In order to do it, we will represent the 
scalar field equation <j> + 3H<j) = — V 1 (fa in the form 

Thus in order to find the total change of the energy density of the scalar field during some 
time one should integrate —3Hfa: 

t 

Ap = A(V + 4> 2 /2) = -3 J Hfadt . (6.10) 

to 

Using this equation, one can find the change of the energy density of the field (j) during the 
time ir/m when the field <p moves from fa to —fa: 



3tt 2 o , 4 ttV3AV 



Ap_ = ^m 2 tf ~ " V ;p mfa{ - ^AVcfi . (6.11) 

In the most interesting case fa « fa, one can neglect the last term in this equation and 
replace fa by fa: 

3vr 2 „ 4 ttVSAV n 
Ap. = -^m 2 ^ - m<p 2 . (6.12) 

Thus, if the initial kinetic energy of the field is equal to zero at the beginning of the 
oscillation at <j) = fa, at the moment when the field <fi will reach the point —fa its kinetic 
energy will be positive, 

— = Ap_ = — m fa mfa . (6.13) 



2 
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Note that for AV < m?(f)Q the last term is much smaller than the first one, so one 
finds, in the first approximation, that the field <j> coming to the point —<j)\ acquires kinetic 
energy 



2 ' 16 

and velocity 



4> 2 , 3vr 2 n 4 3vr 2 ^ 

%r- = Ap_ « -TTrm 2 ^ = — ^ Vb « Vb , (6.14) 



^m0 2 . (6.15) 



This velocity continues to grow during subsequent oscillations and eventually the scalar 
field <j) and the scale factor a blow up, as shown in Fig. |2| 

So far we have studied an expanding universe that stops its expansion and collapses. 
But what if it was collapsing at the beginning of the oscillation? Suppose the scalar field 
was moving very slowly until it reached the point <j)\. Then it started falling down, just as 
in the case considered above. However, this time we will assume that the universe was not 
expanding but collapsing. This corresponds to the choice a = at the beginning 

of the process. 

In this case the universe will continue collapsing with ever growing speed. The evolution 
of the field <f> can be studied by the same methods as the ones used above. The main 
difference will be that the field <j> passing through the point (fr = —<fri ~ (fro W1 U have kinetic 
energy 



_ = Ap + = —m O + m0 o • (6.16) 

The kinetic energy of the field at = — 4>o differs from that at (ft = —(pi by AV^. 
However, for AV -C m 2 (pQ this difference is much smaller than each of the terms in Eqs. 
(6.12), ( 6.16) ). Thus these two equations with the above-mentioned accuracy give the kinetic 



energy of the field not only at <j> = — <j)\ but also at 



90- 



This discussion, as well as the difference between Ap_ and Ap + , will play an important 
role in our investigation of the cyclic universe scenario . As we will see, the cyclic regime 
is possible only if the field (ft, after bouncing from the singularity, approaches the point — (fro 
with energy density greater than Ap + , which in its turn is greater than Ap_, which is the 
energy of this field at the point — (fro on its way towards the singularity. Thus one needs 
this field to bounce from the singularity with an increased energy, and one should check 
that the possible source of this additional energy does not create problems for the scenario. 
In fact, we will see that with an account taken of particle production, the required energy 
increase can be much greater than the difference between Ap + and Ap_. 
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7. Other models with V(4>) < 



Until now we have studied only one simple model with a quadratic potential. However, 
many features of models with negative potentials are model- independent. Consider, for 
example, the model with the "inverted potential" V(<fi) = Vq — m ^ with Vq > 0. This 
is the simplest example of a potential unbounded from below. The evolution of the scalar 
field and scale factor in this model is shown in Fig. 10. As we see, in the beginning the 
universe experiences a stage of inflation when the scalar field slowly rolls from the top of 
the effective potential. (We considered a model with Vq 3> m 2 .) Later on, inflation ends 
and the speed of the field increases. If one neglects the effects of the expansion of the 
universe, at large 6 one has ft = 2(V - V(4>)). Therefore 



a 
a 



1 



{V{cf>) - </> 2 ) = V{<j>) 



-Vn 



(7.1) 



3 V y ' ' VT ' 3 

At large <fi the universe starts moving with ever growing negative acceleration. If one takes 
into account the expansion of the universe, 4> 2 becomes even smaller, and the deceleration 
is even greater. As a result, the expansion slows down and the universe starts contracting. 
At this stage the "friction term" 3H(p in the equation of motion of the scalar field becomes 
negative, which causes the field (ft to grow and leads to a rapid collapse of the universe. 



4> 



a 



Figure 10: Evolution of the scalar field and scale factor in the model V(<j>) = Vq 
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Another example is the standard potential used for the description of spontaneous 
symmetry breaking, with the addition of a negative cosmological constant Vq < 0: 

V{4>) = \{4> 2 - v 2 ) 2 + V = -\m 2 4> 2 + g^ 4 + \m\ 2 + V . (7.2) 

Here m 2 = Xv 2 and the point (ft = v corresponds to the minimum of V{4>) with symmetry 
breaking. The potential V(<f>) becomes equal to Vq < in the minimum of V{4>) at (f> = v. 



As we see in Fig. 11, the scalar field in this case experiences a stage of oscillations near the 
minimum of the effective potential with V(cp) = Vq < 0, but then it jumps off the minimum 
and blows up because of the "negative friction" in the collapsing universe. For most model 
parameters and initial conditions, if the field originally moves towards the minimum with 
cj) = +v it will blow up in the direction <f> — > — oo and vice versa. The reason is that at the 
initial stages of the development of the instability the field 4> is most efficiently accelerated 



- 25 - 



by the negative friction if for a while it moves in a relatively flat direction, i.e. from one 
minimum to another, instead of directly moving upwards ||. 



4> 




a 




Figure 11: Evolution of the scalar field and the scale factor in the model V{<j>) 
with Vq < 0. 



■V , 



When the field accelerates enough it enters the regime 4> 2 V(4>) and continues 
growing with a speed practically independent of V(<p): <ft ~ ln^, 4> ~ see Eq. ( 2.1S| ). 
So for all potentials V((f>) growing at large cf> no faster than some power of 4> one has 4> 2 /2 
growing much faster than V((j)) (a power law singularity versus a logarithmic singularity). 
This means that one can indeed neglect V((f>) in the investigation of the singularity, virtually 
independently of the choice of the potential. Thus we see that from the point of view of the 
singular behavior of the field (f>(t) and the scale factor, potentials having a global minimum 
with V(4>) < are as dangerous as potentials unbounded from below. 

Since a small modification of the potential (shifting the minimum of V{4>) towards 
V{4>) < 0) may lead to a change of regime from expansion to contraction, one may wonder 
whether some other modification of V(4>) can switch the regime of contraction back to 
expansion? The answer follows from the equation H = — ^(p + p). This equation implies 
that H < because p + p > in accordance with the null energy condition. This means, 
in particular, that if the universe switches from expansion to contraction, it cannot later 
return to the regime of expansion. The only possible exception would be if the universe 
were to pass through a stage of super-Planckian density in which the Einstein equations 
were invalid. 

Even though many properties of the theories with negative potentials are model- 
independent, the topology of their phase portraits depends on the choice of the potential 
V(4>). For example, the hypersurface representing a flat universe in the theory V(cp) = 
Vq — — V is given by the constraint equation 



6H 2 + m 2 6 2 



2V n 



(7.3) 



This equation describes a hyperboloid just like the flat universe hypersurface of the theory 
V{<p) = Vq + — g . In this case, however, the axis of the hyperboloid is in the <j> direction 
rather than the H direction. Moreover, the hyperboloid for this model has two sheets 
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for Vo < and one sheet for Vo > 0, which is the reverse of the situation for V(4>) = 

2 J, 2 

Vo + — 2 . The different orientation of the hyperboloid means, for example, that for the 
theory unbounded from below all trajectories end in a Big Crunch singularity, regardless 
of the signs of Vo and k. 

8. Approach to the Singularity, Quantum Corrections, and Particle Pro- 
duction 

Talking about the dynamics of the cosmological scalar field, until now we have remained 
in the realm of classical physics. We ignored possible quantum effects, and in particular 
the effects of particle production. These effects may lead to some important qualitative 
changes of the phase portraits, however, especially near the singularity. 

First of all, near the singularity one may need to take into account quantum corrections 
to the effective action of general relativity. Even ignoring possible effects related to brane 
cosmology or M-theory, one may need to add to the effective action terms proportional to 
R 2 , R^W", etc. 

An important example of such a theory is given by a combination of scalar field theory 
and the Starobinsky model, where the effective Lagrangian has additional terms ~ R? |^4j . 
Whereas this addition is not very significant at low energies, it completely changes the 
behavior of the theory near the singularity. 

For example, in the absence of this term the generic regime for a scalar field approaching 
the singularity is <j> 2 /2 3> V((p), which corresponds to the equation of state p = p. This 
regime was recently discussed in ]29| in the context of string cosmology. As we have seen, 
in this case a ~ t 1//3 , </> ~ lni. 

However, if one adds the term R 2 , the most general regime for theories where the 
potential is not too steep becomes quite different: a ~ i 1 / 2 , 4> ~ i -1 / 2 p4|]. 

It is even more important to consider the effects of particle production. If one ignores 
quantum effects, one typically finds the curvature R ~ t~ 2 in a collapsing universe. Scalar 
particles minimally coupled to gravity, as well as gravitons and helicity 1/2 gravitinos 
|30t| , are not conformally invariant; their frequencies thus experience rapid nonadiabatic 
changes induced by the changing curvature. These changes lead to particle production 
due to nonadiabaticity with typical momenta k 2 ~ R ~ t~ 2 . The total energy-momentum 
tensor of such particles produced at a time t after (or before) the singularity is T^ v ~ 
0(k 4 ) ~ i? 2 ~ r 4 

matter or radiation density of the universe p ~ t~ 2 , one finds that the density of created 
particles produced at the Planck time t ~ 1 is of the same order as the total energy density 
in the universe. 

The main point of this discussion is that particle production near a cosmological singu- 
larity can be extremely efficient. Generically one expects that when the universe emerges 
from or approaches a singularity and its density is close to the Planck density, the density 
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31, 



Comparing the density of produced particles with the classical 



of produced particles should be comparable to the total energy density of the universe. 

This is a pretty general conclusion. For example, in brane cosmology a similar effect 
of particle production may occur even though R = in 4D. Indeed, the change of distance 
between branes leads to a nonadiabatic change of the spectrum of Kaluza-Klein modes and 
thus to particle production; one may call it a time-dependent Casimir effect. Note that 
this effect exists even in theories with unbroken supersymmetry |33| |. 

This observation has many implications. In particular, one can no longer expect that 
matter (or a scalar field) has the equation of state p = p near the singularity. Even if 
the universe around the Planck time was dominated by matter with p = p, the creation 
of particles would immediately change the situation. And even if the density of created 
particles initially was somewhat smaller than the energy density of matter with p = p, 
this situation would rapidly change. The density of the component of matter with p = p 
decreases as a~ 6 , whereas the energy density of radiation and nonrelativistic particles 
decrease as a -4 and a -3 respectively. Therefore the energy density of such particles soon 
becomes greater than the energy density of the matter component with p = p. Once this 
happens the scalar field immediately freezes. It loses its initial kinetic energy and begins 
moving very slowly. As we already discussed, this provides perfect initial conditions for 
inflation. This result also has important implications for the cyclic universe scenario [10]. 



9. Cyclic universe 
9.1 The basic scenario 

Until now we have studied the evolution of the universe and classified new possibilities that 
appear in scalar theories with negative potentials. This problem is very interesting. Its 
investigation has already brought us to an important realization: We cannot live in anti 
de Sitter space dominated by a negative cosmological constant, not because the negative 
cosmological constant is forbidden, but because a universe dominated by negative vacuum 
energy cannot appear after a long stage of inflation |6|, 0, Another interesting realization 
is that the available observational data can tell us nothing about the future of the universe: 
we may live in a stage of a nearly constant de Sitter-like inflationary acceleration, but it 
may end with a global collapse |3j, [35], [!(], |?], || . 

A common feature of cosmological evolution in models with negative potentials is that 
it begins in a singularity and ends in a singularity, even if the universe is not closed. This 
was not the case for the theories with V(4>) > 0, where the universe may continue expanding 
forever and never end in a singularity even if it is closed. 

This naturally brought back old speculations about the oscillating, or cyclic, evolution 



of the universe, see e.g. [16] - [22], [p. Oft . The universe may be created in a singularity, then 
collapse and re-emerge again. 

There is a certain intellectual attractiveness in this idea. However, during the last 20 
years this idea has lost some of its initial appeal. Indeed, if there was a stage of inflation 
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after the singularity, then the initial conditions producing our universe are nearly irrelevant 
for the investigation of the formation of large-scale structure in the observable part of the 
universe. Moreover, inflation in many of its simplest versions is eternal 37, [38|| . This fact 
may not solve the singularity problem [^] , but it puts the origin of our part of the universe 
indefinitely far away in the past 



Recently Steinhardt and Turok proposed a version of inflationary theory where the 
stage of inflation occurs after formation of the large scale structure of the universe and 
perturbations responsible for the formation of the structure of the universe are produced 
before the singularity, during the previous cycle of the universe evolution flO[| . In this 
scenario inflation does not protect us from all uncertainties associated with the physical 
processes occurring around the big bang. On the contrary, in order to describe our universe 
in this scenario one must know exactly what happens with small perturbations of the metric 
when they pass through the singularity. 



The cyclic scenario [1C] is a modified version of the ekpyrotic scenario pl| . It is based 
on the idea that we live on one of two branes whose separation can be parametrized by 
a scalar field <p. It is assumed that one can describe the brane interaction by an effective 
4D theory with the effective potential V((f)) having a minimum at V{4>) < 0. In the 
original version of the ekpyrotic scenario it was assumed that V{4>) is always negative, but 
it vanishes at <p = and at (f) — > oo. It was claimed that one of the main advantages of 
the ekpyrotic scenario was the absence of a cosmological singularity and the possibility 
to solve the major cosmological problems without the help of inflation, which was called 
"superluminal expansion." 

However, later it was found that it is difficult to solve the cosmological problems in 
the ekpyrotic scenario without using inflation [12]. Moreover, perturbations of the field <f> 
that could be responsible for large scale structure formation in this scenario are generated 
due to tachyonic instability [41| at the time when V(4>) was supposed to be smaller than 
— 10 -50 in Planck units. Therefore it is difficult to avoid inflation in this model: Even 
a miniscule positive contribution to V(4>) of the order 10~ 50 would lead to a stage of 
exponential expansion of the universe at large (j> [12|. Also, in 42] it was shown that in 
the context of the effective 4D theory used in [11] the universe can only collapse. This 
means that the ekpyrotic scenario suffers from the cosmological singularity problem. This 
problem has been analysed in 13], but so far it remains unresolved. 

In the cyclic universe scenario the authors assume, in accordance with the suggestion 
of Ref. |i~2|| , that V(4>) is positive at large 4>, and therefore the universe experiences a stage 
of inflation. This stage provides the solution to the major cosmological problems. However, 
it is assumed that this is an extremely low-scale inflation associated with the present stage 
of acceleration of the universe in a state with V((p) ~ 10~ 120 . Inflationary perturbations 
produced at this stage have wavelengths comparable to the present size of the horizon, so 
they cannot be responsible for galaxy formation. 

Therefore it is assumed that the desired perturbations of the scalar field are produced 
after inflation, by the same tachyonic mechanism as in the ekpyrotic scenario [fil], [ll], |l2|] . 



- 29 - 



V 




4> 



Figure 12: Scalar field potential in the cyclic scenario. The minimum of the potential may occur 
at any value of 4>, in this section for simplicity we will assume that it occurs at (\> = 0; we will 
consider a more general situation later. 



The effective potential of the scalar field in the cyclic scenario has the shape shown in Fig. 
12. Inflation occurs at large (p. Once the field rolls down to the region where V(4>) < 0, the 
universe begins to collapse. At that time perturbations of the scalar field are generated. 
The speed of the field in a collapsing universe grows. It reaches the plateau at (j) — > — oo 
where, according to jl(|, the potential vanishes. The universe enters the regime where its 
energy density is dominated by the kinetic energy of the scalar field, and it evolves towards 
the singularity in accordance with Eqs. (2.17), ( |2.18|) . 

Usually, this would be considered the end of the evolution of the universe. However, 
in the cyclic scenario it is assumed that the universe goes through the singularity and 
re-appears again. When it appears, in the first approximation it looks exactly as it was 
before, and the scalar field moves back exactly by the same trajectory by which it reached 
the singularity |]1 



This is not a desirable cyclic regime. Therefore it is assumed in [ 10 1 that the value of 
kinetic energy of the field <p increases after the bounce from the singularity. This increase is 
supposed to appear as a result of particle production at the moment of the brane collision 
(even though one could argue that usually particle production leads to an opposite effect). 
If the increase of the kinetic energy is large enough, the field <j> rapidly rolls over the 
minimum of V(4>) in a state with a positive total energy density, and continues its motion 
at <j) > 0. The kinetic energy of the field decreases faster than the energy of matter produced 
at the singularity. At some moment the energy of matter begins to dominate. Eventually 
(a few billion years after the big bang) galaxies form. Then the energy density of ordinary 
matter becomes smaller than V(4>) and the present stage of inflation (acceleration of the 
universe) starts again. 

As we see, this version of the ekpyrotic scenario is not an alternative to inflation 
anymore. Rather it is a very specific version of inflationary theory. The major cosmological 
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problems are supposed to be solved due to exponential expansion in a vacuum-like state, 
even though the mechanism of production of density perturbations in this scenario is non- 
standard. Let us remember that Guth's first paper on inflation @ was greeted with so 
much enthusiasm precisely because it proposed a solution to the homogeneity, isotropy, 
flatness and horizon problems, even though it didn't address the formation of large scale 
structure. The Starobinsky model that was proposed a year earlier p] could account for 
large scale structure and the observed CMB anisotropy [43], but it did not attract as much 
attention because it did not emphasize the possibility of solving these initial condition 
problems. 

In fact, the stage of acceleration of the universe in the cyclic model is eternal inflation. 
Indeed, the main criterion for the process of self-reproduction of the inflationary universe to 
occur is that the amplitude of inflationary perturbations 5<ft ~ H ~ \/V should be greater 
than the change A(j) of the classical value of the field (j) during the time H : Acj) ~ V' /V 
p7| , j3§| , [iOfl . For the potential V{4>) used in the cyclic model one has = const in the 
limit 4> — > oo, whereas A(p —* in this limit. Thus the universe at large 4> enters the 
stage of eternal self-reproduction, quite independently of the possibility to go through the 
singularity and re-appear again. In other words, the universe in the cyclic scenario is not 
just a chain of eternal repetition, but a growing self-reproducing inflationary fractal of the 



type discussed in |37|, [3§|, |40(] . 



One may wonder, however, whether this version of inflationary theory is good enough 
to solve all major cosmological problems. Indeed, inflation in this scenario may occur only 
at a density 120 orders of magnitude smaller than the Planck density. If, for example, one 
considers a closed universe filled with matter and a scalar field with the potential used 
in the cyclic model, it will typically collapse within the Planck time t ~ 1, so it will not 
survive until the beginning of inflation in this model at t ~ 10 60 . For consistency of this 
scenario, the overall size of the universe at the Planck time must be greater than / ~ 10 30 
in Planck units, which constitutes the usual flatness problem. The total entropy of a hot 
universe that may survive until the beginning of inflation at V ~ 10~ 120 should be greater 
than 10 90 , which is the entropy problem [25|. An estimate of the probability of quantum 



creation of such a universe "from nothing" gives P ~ e I 5 ' ~ exp (—^rr) ~ e 120 44]. 



V 

There are some other unsolved problems related to this theory, such as the origin of 
the potential V{<f)) ll^] and the 5D description of the process of brane motion and collision 



|42| , 45 1 . In particular, the cyclic scenario assumes that the distance between the branes is 
not stabilized. Thus one would need to find some other mechanism that would ensure that 
the effective gravitational constant, as well as other parameters depending on the field (ft 
(i.e. on the brane separation), does not change in time too fast. This is one of the reasons 
why it is usually assumed that the branes in Hofava-Witten theory must be stabilized. 

We will not discuss these problems here. Instead of that, we will concentrate on the 



3 It is remarkable that quantum effects and the mechanism of self-reproduction may work even at the 
present stage when the wavelength of inflationary fluctuations is greater than the size of the observable part 
of the universe and the square of their amplitude is as small as 10~ 120 in Planck units. The reason why it 
may work is that the curvature of the effective potential at large (j> is even much smaller. 
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phenomenological description of possible cycles using the effective 4D description of this 
scenario. This will allow us to find out whether the cyclic regime is indeed a natural feature 
of the scenario proposed in ju]]. 

For the remainder of this section we will analyse this scenario using the tools developed 
in the earlier sections of the paper. In Section |9.2| we will describe the phase portrait of 
the cyclic scenario. In Section EO we will consider the conditions that must be satisfied 



at the bounce in order for the cyclic regime to occur. In Section 9.4 we will analyse 



the motion of the field as it returns from the singularity and show that the conditions 



described in Section |9.3j are difficult to realize self-consistently without invoking super- 



Plan ckian potentials, even in the vicinity of the minimum. Following the authors of 110 



we will consider such super-Planckian potentials in Section 9.5. Aside from the problem of 



applying the effective 4D theory at such high energies, we will find that there are still other 



problems in such realizations of the scenario. In Sections 9.6 and (O we will propose some 



modifications of the cyclic scenario that may resolve some of the problems raised here. 
9.2 Phase portrait of the cyclic universe 

The phase space of the cyclic scenario is the usual 3d space (0, </>, H). If one does not take 
into account matter and radiation, the phase portrait of the scenario forms a 2d surface in 



3d space. It is shown in Fig. 13 without the Poincare mapping. (If one adds radiation, the 
flow of trajectories becomes three dimensional.) The trajectories corresponding to different 
initial values of <fi and <p start at large H, i.e. in the upper part of Fig. [l^. The trajectories 
beginning at large positive <j) reach the (red) separatrix going from the point SI to the point 
A. Its upper part {H > 0) corresponds to inflation. These trajectories follow the separatrix 
towards the throat of the phase portrait at H = 0, and then all of them move towards the 
singularity. The trajectories beginning at large negative <f> fall from the singularity at large 
positive H to the singularity at large negative H without entering the stage of inflation. 

If one flips <fi — > — cj) and H — ► —H, which corresponds to time-reversal, the red 
separatrix connecting points SI and A becomes the green separatrix connecting points R 
and S2. In the lower part of the figure (at negative H) this line corresponds to the stage 
of deflation (exponential contraction of the universe, which is a time-reversal of inflation). 
These two separatrices divide all trajectories into three topologically disconnected parts: 
the trajectories to the right of the green separatrix, the trajectories between the green and 
the red separatrix and the trajectories to the left of the red separatrix. 

One could think that the green separatrix separates inflationary trajectories from the 
trajectories that fall to the singularity without reaching the stage of inflation. However, 
it is not so. As we already discussed in Section ||, the trajectories that reach the stage of 
inflation are at a finite distance to the right away from the green line connecting points R 
and S2 (i.e. at greater values of <fi and (/>). 



The (4>, (f)) projection of the phase portrait for the cyclic scenario is shown in Fig. 14 



also without the Poincare mapping. An interesting feature of the right panel of Fig. 14 



is the apparent absence of any trajectories near the green line (the right separatrix at the 
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Figure 13: The 3d phase portrait for the cyclic scenario. All trajectories (blue lines) begin at 
H > and end in a singularity at H < 0. 



right panel). This might seem surprising because this line is a solution of the equations 
of motion, so there must be other solutions nearby. The reason is that the deflationary 
universe regime described by this line is a strong repulsor, just opposite to the fact that the 
inflationary red line at H > (the right separatrix at the left panel) is a strong attractor. 
As a result, the density of trajectories near the green line at H < is very small; that is 
why they do not show up in Fig. 14. We discussed a similar issue in Section [|. 

As we see, all trajectories beginning at H > end up in the singularity at H — * —oo. 
In the cyclic scenario it is assumed that the universe goes through the singularity and 
re-appears again. When it happens, all trajectories with (ft < 0, (ft < and H < the left 
lower part of the right panel in Fig. 14 suddenly reappear in the right upper corner of the 
left panel of Fig. 14, describing the trajectories starting at (ft < 0, (ft > and H > 0. If one 
ignores particle production at the singularity, the red separatrix on the right panel becomes 
the green line at the left panel (time-reversal). As a result of this flip, the field (ft, which 



previously was running down along the red separatrix towards the singularity in Fig. 13 



eventually returns exactly to the same place at (ft > where it was in the very beginning 
of the process. However, it returns back not at the stage of exponential expansion but at 



the stage of exponential contraction, following the green separatrix in Fig. 13 



Exponential contraction is not a desirable regime. In order to reach the cyclic infla- 
tionary regime, some of the trajectories to the left of the red separatrix after the singularity 
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Figure 14: The 2d phase portrait for the cyclic scenario. All trajectories begin at the bounding 
box of the left panel (H > 0) and end at the bounding box of the right panel (H < 0). 



should jump sufficiently far away to the right of the green separatrix. As we already men- 



tioned, Ref. [10] assumes that this jump may occur due to an increase in the energy of 
the scalar field bouncing back from the singularity. This increase in energy is supposed to 
happen due to particle production. Only if this jump is sufficiently large can these trajec- 
tories reach the red inflationary separatrix going from SI to A. Then inflation begins, the 
field rolls to the minimum of V((f>) again, and everything repeats. 



9.3 Moving towards the minimum of V(4>) 



To study the potential shown in Figure 12 we will assume that near the minimum it can 



be represented as ^-(4> — 4>q). At <fi ^ 4>o we will take it to be flat with V ~ 10 and at 
4> < (fro we will take V = 0. The results of a numerical investigation for more complicated 
potentials are very similar to the ones obtained for this simple model. However, in this 
model one can study everything analytically using the results obtained in Section ^. Indeed, 
we know how the field moves at <j) < — <j)Q, when V(4>) = 0, and we also know how it behaves 
in the quadratic potential, when it moves from —fa to <f>Q. The only thing that we need to 
do is to patch these two regimes together. 

At the initial stage the scalar field moves extremely slowly at <p > 4>o and the universe 
inflates. Once it reaches <j> ph 4>q it falls down, V(4>) becomes negative, and the universe 
begins to contract. To describe this process one can use the theory developed in the first 
part of this paper. The contraction begins at (j) = 4>c fl6.4p. The scalar field reaches <fi = —(po 
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with energy A/j_ given by Eq. ( 6.11 ) 



Subsequently, the field <f> moves towards <f> = — oo and the singularity develops in 
accordance with Eq. ( [2.18 ). To describe this motion one should take to = l/\/3Ap_ in 



(2.18) and replace <^>o by 



In this solution <ft = — <j)Q at t = to- 

Let us use this equation to find the value of the field <j> at the Planck time when the 
energy density becomes 1 in Planck units and one can no longer study this regime within 
the context of general relativity. This happens at t p = l/v3 in Planck units. Therefore the 
scale factor of the universe a ~ t 1 / 3 decreases by a factor ~ (A/9_) 1 / 6 from the beginning 
of the process &t <fi = —<fio until the density becomes O(l). The scalar field cj> at that time 
is given by 

<f> p = -fo + JI lnAp_ . (9.2) 



Setting a(t p ) = 1 we can write our solution as 

0-^ = y|lnV3t, f = ^> a = &*&*, (9.3) 
which in turn implies 

\f2 r- 

</>=-^3-, - <j) p = V6 In a . (9.4) 

One can also represent our results in terms of the conformal time r, where dt = adr. 
In this case t = (2r 3~ 5 / 6 ) 3 / 2 , and 

«^-0 P =y|ln^. (9.5) 

The Planck time t p = l/y/3 corresponds to r p = \/3/2. 

The cyclic scenario requires that the universe bounce back from the singularity and 
the field move back from — oo to <j)Q. Depending on how much kinetic energy the field has 
at this point three regimes are then possible: 

1. ^- < Ap_ at <j) = —0o- This is the regime that would be reached if the bounce 
were perfectly symmetric (in which case 4- = A/)_). The universe starts collapsing 
at (j) < <p c . The field overshoots the point (f> = cpo and moves with ever growing 
speed towards <j> = +oo. There is a small bunch of trajectories such that the scalar 
field evolves very slowly, the equation of state is p = —p, and the universe contracts 
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exponentially. Eventually, however, the kinetic energy of the field (j) dominates and 
the collapse becomes power-law with p = p. 

This regime is represented by the trajectories to the left of the green separatrix in 



the upper part of the left panel in Fig. 14 



2. Ap_ < 4r < Ap + at <fi = — 0o- The universe starts collapsing at 4> > 4>c- The field 
does not have enough energy to reach the point <fi = <pQ, so it returns back to negative 
(f>, the field moves with ever growing speed to <f> = — oo, and a singularity develops. 

This regime is represented by a small bunch of trajectories to the right of the green 
separatrix in the upper part of the left panel in Fig. [l4|. 

3. 4p p> Ap + at <p = — (j)Q. The universe continues expanding and the field (f> becomes 
greater than <f>Q. It continues growing and gradually slows down. As a result, inflation 
begins. Then the field very slowly decreases, falls into the minimum of V(cj)), the 
universe collapses and the field moves to (j) = — oo. This is the regime required by 
the cyclic scenario. 

This regime is represented by trajectories starting sufficiently far from the green 



separatrix, to the right of it in the upper part of the left panel in Fig. 14 



The last of these regimes requires additional explanation. Let us remember how we 
derived the expression for Ap + : We considered the field <f> slowly rolling from = 0o 
during the stage of contraction and found that it arrived at the point <fi = —<j)Q with kinetic 
energy Ap + . If we reverse the time evolution of the universe, we will see the scalar field 
rolling down from 4> = — 4>q and arriving at the point 4> = 4>q with a nearly vanishing speed 
during the stage of expansion. If the initial kinetic energy of the field is greater than Ap + , 
it reaches the point <f> = (f>o with a non-vanishing speed and moves further onto the plateau 
where the energy density of the field 4> becomes constant, and inflation begins. 

As we have seen in Section ||, the difference between Ap + and A/?_ is extremely small: 

bp = Ap + - A/?_ = irV3AVm<pl . (9.6) 

Here AV has the meaning of the height of the effective potential at (f> > 0o> m ° ur case 
AV ~ 10~ 120 . Thus one might expect that it pretty easy to jump from the trajectory with 
energy Ap_ to the desirable trajectory with energy greater than Ap + , as in case 3. 

In reality, however, the required jump in kinetic energy becomes much larger when 
one takes into account quantum effects. As the field eft moves through the minimum from 
— (j)Q to 0o its mass changes from to m and back to again, all within a time 0{m~ l ) 



(half of an oscillation), see Fig. 1£. This non-adiabatic change, 4r? ~ m2 > wm l ea d to the 



production of ^-particles with energy density 0(m 4 ) |27]]. Therefore the field (j) loses an 
amount of energy 0(m ), which makes it less likely to reach <fio while the universe is still 
expanding. 4 Thus in order to realize the cyclic scenario the kinetic energy density of the 

4 Note that the production of ^-particles during this very short time interval appears in addition to the 
process of particle creation near the singularity discussed in Section [| 
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field <j> at the point — 0o must be greater than Ap + by 0(m 4 ), which is much greater than 

Ay. 

One may wonder where the field gets this boost in kinetic energy from. Usually one 
would expect that the field after a bounce can only lose energy due to particle production. 
However, in [ 10 1 it is assumed that it can actually gain energy as a result of particle 



production during the brane collision (i.e. in the singularity). It is not quite clear whether 



this can indeed happen, see e.g. [45| where it is claimed that particles can be created 



during the brane collision only if they have negative energy density. We are not going to 




discuss this issue here. Instead of that, we will follow the assumptions of [1C] and check 
what happens to the scalar field (p if the universe after the bounce contains some matter 
or radiation. 

9.4 A scalar field with a vanishing potential in the presence of radiation 

Let us consider the motion of the field </> from — oo to — </>o in the presence of radiation. 
The Friedmann equation describing this process can be written as follows: 

2 

Here <pi is the velocity of the field at some moment t{ , dj is the scale factor of the universe 
at that moment, and p\ is the density of radiation at that time. This equation reflects the 
fact that the kinetic energy of the field decreases as a -6 and radiation energy decreases as 
a -4 during the expansion of the universe. 5 

It is convenient to write this equation in terms of the conformal time r, where dt = adr: 

(a') 2 = ^ + B, (9.8) 

where a' = ^ = ad, A 2 = <p 2 a 6 j6 and B = 

Taking a(0) = (at the singularity), the solution of this equation is 

a 2 = 2At + Bt 2 . (9.9) 

For definiteness, we will normalize our solution at the time U = t p , when (j) 2 /2 = 1 and 
at = 1. Then A 2 = 1/3, B = p£/3, and 

2= 2 ^2 

^3 3 V7 



Then, using equation <j)' = ^a?/a 2 = V6A/a 2 , one finds 



| In ^ 2T - r + C r = ^ln °-^ + C r . (9.11) 

2 VSCl + ^r) l + ^r 



Here we are considering processes at sub-Planckian energies where the usual Friedmann cosmology is 
supposed to be valid. 
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Here <ft p is the value of the scalar field at the time when (ft jl = 1 after the bounce. The 
constant of integration C r is supposed to vanish in the absence of radiation, i.e. for p p = 0. 
In this case (ft p = <p p , and our solution ( |9.11| ) coincides with the solution presented in Eq. 
Q9.5| ). This means that in the absence of radiation the field (ft elastically bounces from the 
singularity, in accordance with [13]. 



One can find the constant C r for any given p p from the condition that (ft = (ft p at 
b 2 /2 = 1 and a = 1. In particular, for p p <C 1 one has C r ~ 

Eq. fl9.ll ) implies that at p p r > 2\/3 the field stops moving. Therefore we will assume 



that PpT <C 1 at (ft < —(f>o- This leads to a strong constraint on p p : 

P ; & (Ap^ 3 . (9.12) 

If one takes, for definiteness, <po ~ 0.1M P , m 2 ^ ~ 10~ 20 , as in the original version of the 
ekpyrotic scenario [11], one finds that the cyclic scenario with these parameters cannot 
work unless the energy density of radiation at the Planck time is less than 10 -6 in Planck 
units. In general the density of gravitationally produced particles is ~ H 4 , which is O(l) at 
the Planck time, so it is not clear how particle production could be so strongly suppressed. 

Suppose, however, that for whatever reason one can indeed have p r p <C (Ap^.) 1 / 3 . In 



this case p r r <C 1 and Eq. (9.11) can be represented in the following form: 



V61na-^(r-^). (9.13) 



With our normalization of a one has 



\fi = a" 6 , (9-14) 



As we already discussed, if we want the field to move to (ft > (fro during the stage of 
expansion of the universe, its kinetic energy (ft 2 /2 must be greater than Ap + at (ft = —(fto- 
If we assume that the field has sub-Planckian energy as it moves through the minimum, 
i.e. that Ap + <C 1, then 

4> p >- ( j )0 + — In A P+ + OL^ (Ap+y 1 / 3 . (9.15) 



Comparison with Eq. (9.2) gives the following condition: 



^_, p> ^ (A , + r V3 + J_ ln ^±. (, 16) 

In general, it could happen that after bouncing from the singularity the field (ft appears at 
the Planck density at (ft p 7^ (ft p , so that (ft p — (ft p = 0(p p ) |Io|] . However, our investigation 
shows that the cyclic scenario with Ap + <C 1 could work only if (ft p — <ft p 3> p p . 

This means that the cyclic scenario can work only if a very small amount of radiation 
can produce a major change in the state of the field (ft at the Planck time: <ft p — (ft p ;> 



- 38 - 



Pp (Ap + ) -1 / 3 . Second, the amount of radiation at the Planck time must be very small, 
Pp (A/9+) 1 / 3 . This may be a real problem if, as we expect, quantum effects at Planckian 
densities create particles with density p r p = 0(1). 

These problems are less serious in models with Ap+ > 1, i.e. if the field <p acquires 
super-Planckian energy even before it reaches the plateau at <j> < — 0o- Such models are 
suspect because the usual 4D approach based on general relativity becomes unreliable at 
super-Planckian densities. It appears that such models are necessary for the cyclic model, 
however, and in at least one of their papers the authors of [|n]] invoke such a model. We 
therefore consider such potentials here. 



9.5 Super-Planckian potentials for the cyclic scenario 

Let us now consider a potential proposed by the authors of the cyclic scenario |T( 

V(cP) = V (1- e- c *) F(0) . (9.17) 



-l4> 



In the particular example studied in the last paper of Ref. |T(| one has F((f>) = e~ 
Vq = 10~ 120 , c = 10, 7 ~ 1/8. This potential is shown in Fig. |l^. This potential has 
the same structure as the potential shown in the Fig. 12, but the scales and the position 



of the minimum are determined by the parameters given in [10]. At <f> = this potential 
vanishes. It approaches its asymptotic value Vo = 10 -120 at <j> ^ 1. Inflation in this scenario 
is possible at 4> ^ 1. At ;> 15 one has V ?J / 2 £ V and the universe enters the process 



of eternal inflation [37, 38]. The potential has a minimum at cf) ~ —36; the value of the 



potential in this minimum is V m { n ~ —3. 
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Figure 15: An example of cyclic scenario potential used in Ref. |i0(| . 
Let us try to understand the origin of the parameters c = 10, 7 ~ 1/8 used in |l0| 



According to [14], the amplitude of density perturbations in this scenario in the limit c S> 1 
can be estimated as 



6p 



10-5 yz^ ^ 



(9.18) 
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where Vj is approximately equal to the value of the potential in its minimum V m i n and 
£ is the efficiency with which radiation is produced at the singularity; it is assumed that 
£ <C 1- This suggests that in order to be consistent with observational data ~ 10~ 4 ) 
one should have — Vj 3> 1. This means one must rely on calculations using the equations 
of general relativity at | V(4>)\ » 1. 

The authors of |14| have warned the readers that their results are very preliminary and 
many authors do not agree with their derivation of the amplitude of density perturbations 
1 15]. Therefore it may happen that the correct equation for perturbations in the cyclic 
scenario as well as the expression for V(4>) will be quite different. Here we will simply try 



to understand the values of the parameters used in [1C] and check the consequences of the 
potential they suggested. 

The spectrum of density perturbations obtained in |l4[] is not blue, as in |11], but red, 
like in the pyrotechnic scenario and in the simplest versions of chaotic inflation. The 
spectral index is n ~ 1 — 4/c 2 . Observational data suggest that n = 0.93 ± 0.1, which 
implies that c ;> 5. If one takes c > 5, and Vj > 1, one finds that the curvature of the 
effective potential in its minimum becomes much greater than 1. 

Once one takes V ~ —3 in the minimum of the potential with c = 10 [10], the parameter 
7 can be determined numerically: 7 = 0.1226. It would be hard to provide explanation 
of the numerical value of this parameter. Meanwhile if one takes 7 = 1/8 = 0.125, one 
finds V ~ —3 x 10~ 3 in the minimum of the potential. This would reduce ^ by a factor 
of 30. Thus, in order to have density perturbations with a correct magnitude one should 
fine-tune the value of 7 = 0.1226 with accuracy better than 1%. 




Figure 16: Effective mass squared m 2 
V((f>) in the cyclic scenario. 



V" of the scalar field in the vicinity of the minimum of 



Fig. 16 shows the effective mass of the field (f>. As we see, 



|m 2 | = \V"\ ^ 1 in the 
vicinity of the minimum of the effective potential. A numerical investigation of the motion 
of the field moving from (j) > in a theory with this potential shows that its kinetic energy 
at the moment when (p reaches the minimum of the effective potential is O(10 2 ). When the 
field approaches (ft ~ —39, where the effective potential becomes flat, the kinetic energy of 
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the field (f> becomes ~ 10 6 , i.e. a million times greater than the Planck density. 

Even if we continue to trust our calculations in such a regime, there are still problems. 
First of all, there is a distance A(j) > 30 from the point <fi ~ —30 where the field emerges 
from the deep minimum of its effective potential to the region <f> > 1, where inflation in 
this theory may begin. Let us assume that the kinetic energy of the field is smaller than 
the Planck energy at ~ —30, since otherwise we just cannot trust our analysis at all. 
This assumption is in accordance with |J10|]. Indeed, according to the estimates made in 
|1C], 4> p — 4> p In ^^n) < l m | < 1- I n this model 4> p —34, so indeed one expects 
(b n < —33. 



As we discussed in Section |9.4| , we expect that gravitational particle production will 
create particles with density O(l) at the Planck time. Independently of gravitational 
production, however, there should be production of (f> particles with density 0{m ) due to 
the nonadiabatic change of the effective mass of the field moving from (j) = — 39 to <j> > —32, 
see Fig. 16. In this model 0(m 4 ) ^ 0(1). Thus, when the field reaches the relatively flat 
region at (f> > —32, its motion produces ultra-relativistic particles (j) with super-Planckian 
energy density. These particles, just like usual radiation, immediately freeze the motion of 
the field <j). One can show that in this scenario the field 4> can reach the inflationary regime 
at <p > (which is necessary for the consistency of the cyclic scenario) only if at 4> ~ — 32 
(i.e. at the flat part of the potential) the kinetic energy density of the field eft is 12 orders 
of magnitude greater than the (Planckian) energy density of the produced particles. The 
effective 4D description in terms of the scalar field </> and its effective potential V(4>) is 
inapplicable for the description of such processes. 

This problem is not unresolvable. For example, one may consider effects related to 
non-relativistic particles produced at the singularity. These particles contribute to the 
equation of motion for the field (f> by effectively increasing its potential energy density JlO]] . 
They may push the field towards positive values of the field <f> despite the effects described 
above. However, this would add an additional epicycle to a scenario that is already quite 
speculative. Indeed, one would need to produce a sufficiently large number of such particles 
and make sure that massive particles decouple from the scalar field at the present epoch. 
The last condition is necessary to avoid a rapid change of the coupling constants related 
to the brane separation described by the field 4>. 

One may try to improve the situation by altering the shape of the potential. First of 
all, the original argument of |l(J was that the function F(4>) appears because at small values 
of the string coupling g s nonperturbative effects should be suppressed by a factor e~ 1 ^ 9s or 
e _1 /sf ; or perhaps by e _87r2//9 = . In the case of type IIA (or heterotic) string theory in d=10 
the string coupling is g s = |l3|]. Thus one could expect the suppression function to be 
one of the three proposed types: F(<p) ~ e~ e 7 , F{(f)) ~ e~ e 7 , or F((fi) ~ e~ 8w e 7 , 
with 7 = 1 rather than with 7 = 0.1226. 

It is possible to have V m \ n = —3, as in [ 10 1 , for 7 = 1, but only if one takes c = 81.56. 
The value of c must be fine-tuned: a change in c of 1% results in a change of V m \ n by two 
orders of magnitude. In accordance with |[4|], this would lead to an order of magnitude 
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change in the amplitude of density perturbations. 

With these parameters, however, the curvature of the effective potential in its minimum 
becomes two orders of magnitude greater than the Planck mass squared, so all calculations 
in such models in the context of the effective 4D theory are unreliable. In potentials with 
F((p) ~ e~ e 27<#> or F{4>) ~ e~ Sn2e 7< ^ the curvature in the minimum with 1^(^)1 ^ O(l) 
becomes much greater still. 



9.6 Bicycling scenario 

Various modifications to the cyclic scenario are possible. For example, instead of the 



asymmetric potential shown in Figs. 12 and 15, one may consider a symmetric potential, 



as in Fie. 17 



V 




4> 



Figure 17: Symmetric scalar field potential in the new cyclic scenario. At large values of \(j>\ one 
has V(cf>) « Vb ~ 10~ 120 and there is a minimum at <f> = 0. 

In the beginning, the scalar field is large and positive and it slowly moves towards the 
minimum. When it falls to the minimum the universe begins to contract and the field is 
rapidly accelerated towards the singularity at (ft = — oo. As we already mentioned, the 
structure of the singularity is not sensitive to the existence of the potential, especially if 
it is as small as Vq ~ 10~ 120 . Suppose in the vicinity of its minimum the potential is 



approximately quadratic, V{4>) ~ m 2 {4> 2 
Eq. (EO) the kinetic energy of the field 



<fio)/2. If 0o < 1 and mCl, then according to 



reaches the Planck value at 



+ 



- lnAp_ 



— In (mcfo) 



(9.19) 



For definiteness, suppose that m ~ y/Vo ~ 10~ 60 , and ipo = 0(1). Then we would not even 
know that such a minimum exists (the field would not move there) until the energy density 
of matter dropped below its present density 10 -120 . In this case the kinetic energy of the 
field moving towards (p = — oo would reach the Planck value at (j) p ~ — 112. At that time 



the scale factor of the universe would decrease by a factor of Ap 



1/6 



10 



-20 
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Now let us assume, as in [10], that the field (f> bounces from the singularity and moves 
back. Its energy density drops down to the Planck energy density at <p p ~ <j) p ~ 10 2 . 
During its subsequent evolution the kinetic energy of the field rapidly drops down because 
of radiation. Even if the density of radiation at the time when <j) = cp p were as small as 
10~ 39 , it would eventually begin to dominate because its relative contribution grows as a 2 , 
i.e. up to 10 40 times before it reaches — 4>q. 

Therefore the field 4> freezes at large negative <fi. At this stage the energy density is 
dominated by particles produced near the singularity and density perturbations prepared 
during the previous cycle lead to structure formation. Then the universe cools down while 
the field is still large and negative and the late-time stage of inflation begins. During this 
stage the field slowly slides towards the minimum of the effective potential and then rolls 
towards the singularity at 4> ~^ 00 ■ When it bounces from the singularity, a new stage of 
inflation begins. The universe in this scenario enters a cyclic regime with twice as many 
cycles as in the original cyclic scenario of Ref . |n| . One may call it the bicycling scenario. 

An advantage of this scenario is that it may work even if a lot of radiation is produced 
at the singularity and the field <j) rapidly loses its kinetic energy. However, if in order to 
have density perturbations of a sufficiently large magnitude one needs to have a potential 
with a super-Planckian depth V(</>) < —1, as in @, H, then this scenario has the same 
problem as the scenario considered in the previous section. The kinetic energy of the field 
(f> becomes greater than the Planck density as soon as it rolls to the minimum of V((f>). It 
becomes even much greater when the field rolls out of the minimum, and the 4D description 
fails. 



9.7 Cycles with inflationary density perturbations 

As we see, one of the main difficulties of the cyclic scenario is related to the non-inflationary 
mechanism of generation of density perturbations. It requires a very specific and fine-tuned 
potential, see [12] and discussion above. According to jl4|] , this potential must have a super- 



Planckian depth, so one cannot study the corresponding processes by traditional methods. 
Moreover, the very existence of this mechanism of generation of density perturbations 



remains controversial [15|. 



This problem can be avoided if we consider a potential that grows at large \4>\, such 
as the one shown in Fig. [l^. The field begins to move from large positive (f>, falls to the 
minimum of V((ft), and moves with ever growing speed to —cj). If, for example, the potential 
grows like (fi n at a sufficiently large negative (ft, it does not affect the motion of the field <fi 
towards the singularity. However, when the field <j) bounces back, it immediately loses its 
velocity due to the impact of radiation created at the singularity. Therefore it slows down 
and enters a stage of inflation. At this stage all good and bad memory about the previous 
life of the universe and processes at the singularity are erased and new density fluctuations 
are produced. All particles produced at the singularity become diluted, but new ones 



are produced at the end of inflation due to gravitational effects [46] or by the mechanism 



of instant preheating [47, |4q| . These new particles constitute the matter contents of the 
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observable universe. 
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Figure 18: Scalar field potential in the cyclic scenario incorporating a stage of chaotic inflation. 
Inflationary perturbations are generated and the large-scale structure of the universe is produced 
at 6 < 0. 



Gradually the density of ordinary matter decreases, and the energy density of the 
universe becomes determined by V(4>) ~ Vq. The universe enters a stage of low energy 
inflation (quintessence), which may result in a regime of self-reproduction if V(4>) is flat 
enough. In those exponentially large domains of the universe where the field eventually 
falls down to the minimum of V ((/>), it continues rolling to eft = oo, bounces back after the 
singularity, slows down due to radiation, experiences low-energy inflation, and rolls down 
to the minimum of V(4>) again. 

In this model of the oscillating universe one can have large scale structure formation 
due to inflationary perturbations without any need to rely on controversial assumptions 
about the behavior of perturbations passing through the singularity. Also, one no longer 
needs to have potentials with |V(0)l > 1. However, in this model inflationary perturbations 
are generated only every second time after the universe passes the singularity (at eft < 0, 
but not at (f> < 0). The model can be made even better by making the potential rise both 

In this case the stage of high-energy inflation 



at 



oo and at 



-oo, see Fig. 19 



and large-scale structure formation occurs each time after the universe goes through the 
singularity. 

Thus we see that it is possible to propose a scenario describing an oscillating infla- 
tionary universe without making any assumptions about the behavior of non-inflationary 
perturbations near the singularity. Another important advantage of this scenario is that 
inflationary cycles may begin in a universe with initial size as small as 0(1) in units of the 
Planck length, just as in the standard chaotic scenario 0]. Still, in many other respects this 
scenario is almost as complicated as the cyclic scenario of Ref . |l0| . The theory of reheat- 



ing of the universe in this model, just as in [10], is rather unconventional. Gravitational 
particle production, which is the only source of matter in this scenario, may dramatically 
overproduce gravitinos and moduli fields [46, |48|| . To avoid this problem one would need to 



-44- 



V 




Figure 19: Scalar field potential in the cyclic scenario incorporating a stage of chaotic inflation. 
Inflationary perturbations are generated and the large-scale structure of the universe is produced 
both at cj> < and at <fi > 0. 

use the mechanism of instant preheating |47], |48fl . In order to combine the stage of chaotic 
inflation and the stage of low-scale inflation (quintessence) the potential must be rather 
complicated. To avoid this complication one may need to consider two-field models of the 
type of hybrid inflation. 
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Figure 20: The scalar field potential that appears after the step-by-step simplification of the cyclic 
scenario. 

The main problem of this model is that one still must assume that somehow the 
universe can go through the singularity. But now this assumption is no longer required for 
the success of the scenario since the large scale structure of the universe in this scenario 
does not depend on processes near the singularity. This allows us to remove the remaining 
epicycles of this model. Indeed, the main source of all the problems in this model is 
the existence of the minimum of the effective potential with V((f>) < 0. Once one cuts 
this minimum off, the potential becomes extremely simple, see Fig. and all problems 
mentioned above disappear. In particular, one may use the simplest harmonic oscillator 
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potential ^-(fi 2 + Vq with Vo ~ 1CP 120 considered in the beginning of our paper. This 
theory describes an eternally self-reproducing chaotic inflationary universe, as well as the 
late stage of accelerated expansion (inflation) of the universe driven by the vacuum energy 
V >0. 



10. Conclusions 

The main goal of our work was to perform a general investigation of scalar field cosmology in 
theories with negative potentials. We have found that the phase portraits of such theories 
in the 3d space ((f), (ft, H) have different geometry as compared with phase portraits in 
theories with V((ft) > 0. In theories with V((ft) > the phase portraits for flat universes are 
divided into two disconnected parts describing expanding and contracting universes (H > 
and H < 0). Meanwhile in theories with V((ft) < these two parts become connected. The 
trajectories moving towards V((ft) < simultaneously move from the parts of the phase 
portrait with H > towards the parts with H < 0. Once the universe begins to contract, 
it never returns to the stage of expansion until it reaches the singularity. 

This does not mean that theories with negative potentials should be banned from 
consideration. In some cases the scalar field may be trapped in a metastable minimum, or 
it may roll towards V((ft) < extremely slowly. However, it is quite interesting that with 
an account taken of general relativity potentials that have minima at V((ft) < can be as 
dangerous as potentials unbounded from below. 

A general feature of all trajectories bringing the universe towards the singularity is 
that in all theories with power-law potentials the kinetic energy <ft 2 /2 becomes much greater 
than V((ft) near the singularity. This means that the description of the singularity is nearly 
model-independent, at least at the classical level. In particular, the equation of state of 
the universe approaching the singularity typically is p = p. 

However, this conclusion can be altered with an account taken of quantum effects, 
including particle production near the singularity. Typically particle production near the 
singularity is so efficient that it turns off the regime p = p when a contracting universe 
approaches the Planck density. The effects related to particle production are especially 
significant in an expanding universe as they tend to completely eliminate the stage with 
p = p. 

In addition to the general study of cosmology with negative potentials, we performed an 
investigation of a possibility that our universe may experience repeated cycles of inflation 
and contraction pQ ]. For a complete study of this scenario one would need to resolve 
the singularity problem, as well as several other problems discussed in |4^, 15, 45]. 



In addition, as we show in this paper, the parameters of the effective potentials used in 
the cyclic scenario must be fine-tuned with accuracy better than 1%. This scenario, as 
proposed in pCfl , requires investigation of an effective potential V((ft) of a super-Planckian 
depth, |^(</>)| > 1, and of a scalar field with mass greater than the Planck mass. Even 
if all of these problems could be resolved in the context of a more general approach, the 
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existence of a cyclic regime in the model of Ref. |10[ would require additional assumptions. 
We have shown that ultrarelativistic particles produced near the singularity, as well as 
scalar particles created when the field falls down to the minimum of the effective potential, 
tend to halt the motion of the classical field eft, which prevents inflationary cycles from 
occurring. One way to address this problem is to study quantum creation of supermassive 
particles with specific interactions with the scalar field. However, this would add new 
"epicycles" to a scenario that is already very complicated. 

We proposed several modifications to the cyclic scenario of Ref. that could make 
it more realistic and less dependent on the unsolved singularity problem. In particular, 
if one assumes that the potential V((f)) slowly grows at large \4>\ then the universe may 
still enter a regime of eternal oscillations, but the singularity will be separated from the 
stage of large scale structure formation by a stage of chaotic inflation. This scenario allows 



us to combine attractive features of the oscillating universe model |L(| - [2C] and chaotic 
inflation Q. An important advantage of this model is that it does not need to rely on the 
controversial theory of density perturbations passing through the cosmological singularity. 

But even this model remains very complicated. Fortunately, it allows for one final 
simplification that resolves all of its remaining problems. If one removes the minimum 
of the potential at V(<p) < 0, one returns to the usual scenario of chaotic inflation. It 
describes an eternally self-reproducing inflationary universe, as well as the present stage of 
accelerated expansion. 
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